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ABSTRACT 


A  combination  electrostatic  classifier  -  nucleus  counter  was  used  to 
examine  aerosols  from  three  laboratory  sources.  A  method  of  data  reduction 
was  devised,  not  requiring  computational  facilities,  which  appears  adequate 
if  the  aerosol  is  not  excessively  heterogeneous. 


RESUME 


Une  combinaison  de  classif icateur  4lectrostatique  -  compteur  & 
nucl4ation  fut  utilis4e  dans  le  but  d’ examiner  des  aerosols  provenant  de 
trois  laboratoires.  Une  m4thode  de  r4duction  des  donn4es  a  4t4  d4velopp4e 
ne  n4cessitant  pas  1* utilisation  d’un  ordinateur  et  celle-ci  semble 
ad4quate  dans  le  cas  oil  l’a4rosol  n'est  pas  excessivement  h4t4rog4ne. 
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INTRODUCTION 


This  technical  note  describes  the  use  of  two  instruments,  the 
electrostatic  classifier  (TSI  3071)  and  the  condensation  nucleus  counter 
(CNC,  TSI  3020),  in  the  examination  of  aerosol  size  distributions.  The 
classifier  separates  the  aerosol  into  a  series  of  narrow  fractions  on  the 
basis  of  their  electrical  mobilities,  and  the  CNC  acts  as  a  count 
.instrument  ^for  the  successive  fractions  Kim  Ip© r  t;S>’  Q4  O  1  o.  ^  • 
•*->  ‘  ^  V*-  *  °i-'  '  -•  r'lOl  Ccv 

The  classif ier-CNC  combination  was  used  to  study  aerosolsgenerated 
by  three  sources: 


(1)  The  TSI  3076  atomizer.  This  device,  first  described  by  Liu  and 
Lee  (1)  is  a  development  of  the  Collison  atomizer  with  provision  to 
ensure  long  term  stability  of  output.  Three  of  these  are  integral 
parts  of  the  TSI  FTS  *400  system  in  use  at  DREO.  For  certain 
out-of-system  applications,  it  was  necessary  to  characterize  the 
aerosol.  In  this  work  a  solution  of  0.2J  w/w  dioctyl  phthalate  (DOP) 
in  isopropanol  was  used  as  source,  the  isopropanol  evaporating  from 
the  aerial  suspension  to  leave  an  aerosol  of  DOP  droplets. 

(2)  The  Q127  penetrometer  (Air  Techniques  Incorporated,  Baltimore, 
Maryland).  This  instrument  is  used  primarily  for  canister  testing, 
and  is  based  on  a  technique  first  developed  by  Sinclair  and  La  Mer 
(2).  Liquid  DOP  is  atomized  and  the  original  polydisperse  spray  is 
evaporated  and  recondensed  to  form  an  aerosol  of  considerably  narrower 
size  distribution.  The  output  is  claimed  to  be  droplets  of  0.3  pm 
average  diameter,  at  a  concentration  of  100  ±  20  ug/1.  This  is  either 
a  mass  median  or  mass  mean  diameter;  which  of  these  is  not  clear  from 
the  literature  available  on  the  instrument.  In  any  case  the  system 
had  been  in  continuous  use  with  little  maintenance  for  many  years,  and 
this  work  was  a  check  to  see  whether  its  output  still  met 
specification  requirements. 

(3)  The  model  260  BC  Respirator  Fit  Test  System  (Dynatech  Frontier 
Corporation,  Albuquerque,  New  Mexico).  Here,  liquid  corn  oil  is 
atomized,  and  with  no  further  treatment  apart  from  removal  of  large 
droplets  by  an  impactor,  is  taken  up  in  an  air  stream  and  circulated 
through  a  test  chamber.  The  size  distribution  then  is  expected  to  be 
broader  than  for  the  Q127.  The  mass  mean  diameter  is  stated  to  be  in 
the  range  0.5  -  0.7  pm,  and  the  concentration  with  normal  settings, 

25  ±  5  ug/1. 
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It  had  been  proposed  to  locate  the  Dynatech  chamber  in  a  cold  room 
and  operate  it  at  temperatures  ranging  down  to  -20°C.  The  controls, 
aerosol  generation  and  measurement  systems  would  be  of  course  in  an 
adjoining  room  at  ambient  temperature,  with  air  supply,  return  and  sampling 
lines  passing  through  ports  into  the  cold  room  and  thence  to  the  test 
chamber  itself.  Comparative  tests  with  the  chamber  at  ambient  conditions 
and  under  refrigeration  were  carried  out  to  see  the  effect,  if  any,  on  the 
aerosol . 

In  addition  to  this,  a  brief  examination  was  made  of  the  laboratory  . 
air  furnished  by  the  building  HVAC  system.  Results  are  only  preliminary ' 
and  semi-quantitative;  one  observes  the  expected  dependence'on  the  level  of 
general  activity  in  the  building,  in  the  submicron  range,  but  there  is  an 
apparent  increase  in  concentration  as  the  size  approaches  1  pm,  which  is 
unexpected  and  requires  further  investigation. 

Several  instruments  have  been  described  which  measure  aerosol  size 
distributions  (3~5).  They  all  employ  a  mobility  analyzer,  such  a?  the  TSI 
3071  classifier  referred  to  above.  The  count  instrument  first  employed  was 
the  electrical  aerosol  analyzer  (6),  but  the  CNC  has  more  recently  come 
into  use. 

When  voltage  and  air  flows  have  been  established  for  the  classifier, 
it  will  select,  from  the  incoming  aerosol,  particles  In  a  narrow  range  of 
electrical  mobility.  Uni  positively  charged  particles  in  this  mobility 
range  have  a  correspondingly  narrow  range  of  particle  diameter.  Several 
stages  of  calculation,  however,  lie  between  the  raw  count  data  for 
successive  voltage  settings  and  the  information  ultimately  sought,  number 
concentration  in  each  of  the  diameter  intervals  spanning  the  aerosol.  The 
main  complication  is  that  aerosol  particles  exist  in  a  number  of  charge 
states,  so  that  with  any  given  set  of  operating  conditions,  voltages  and 
air  flows,  the  classifier  will  select  not  merely  singly  charged  particles 
of  one  size  and  the  corresponding  mobility,  but  also  a  fraction  of  doubly 
charged  somewhat  larger  particles,  which  possess  the  same  mobility;  and 
likewise  for  triply  charged,  etc.  Computer  programs  have  been  developed  to 
treat  the  data  by  applying  iteratively  a  series  of  corrections.  The  TSI 
company  has  marketed  a  system,  the  differential  mobility  particle  sizer 
(DMPS,  ref  7)  consisting  of  the  classifier  tnc.  together  with  auxiliary 
equipment  and  a  programme  which  executes  a  series  of  roeasureme.,tz,  setting 
the  required  voltages  and  calculating  distributions. 

In  general,  the  broader  the  dispersion  the  more  extensive  the  data 
treatment  needed  to  recover  the  size  distribution.  It  was  found  here  that 
a  series  of  not  very  laborious  hand  calculations,  to  correct  for  doubly  and 
triply  charged  particles,  would  give  a  quite  reasonable  approximation  to 
the  distribution,  higher  order  corrections  being  negligible.  This 
procedure  was  applied  first  to  the  Q127  and  TSI  3076.  Their  outputs  a~e 
reasonably  homogeneous  and  centred  around  a  size  (several  tenths  of  a 
micron)  known  to  be  well  within  the  capacity  of  the  classifier  at 
reasonable  air  flows,  so  that  the  entire  aerosol  should  be  accessible  to 
measurement.  It  was  anticipated  that  the  Dynatech  output  might  not  be 
equally  amenable  to  treatment. 
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EXPERIMENTAL 


AEROSOL  SAMPLING 


Aerosol  was  sampled  from  the  Dynatech  test  chamber  at  atmospheric 
pressure.  With  the  Q127,  the  hose  conducting  test  aerosol  to  the  canister 
chuck  was  detached,  and  connected  to  the  sample  (polydisperse)  inlet  of  the 
classifier.  The  Q127  is  an  underpressure  system,  air  being  drawn  by  vacuum 
pump,  and  it  appeared  probable  from  the  instrument  manual  that  pressure  at 
this  point  is  near  atmospheric.  Sample  flow,  as  will  appear,  is  only  a 
small  fraction  of  the  main  (sheath  air)  flow,  which  is  drawn  directly  from 
the  ambient  through  a  low  pressure  drop  filter,  so  that  the  classifier  is 
operating  at  atmospheric  pressure  in  this  case  also. 

The  TSI  atomizer  was  sampled  through  the  large  bore  tygon  tube  which 
in  normal  operation  of  the  FTS  HOC  conducts  aerosol  to  the  test  chuck. 

This  line  passes  three  large  capacity  filters  which  introduce  dilution  air, 
and  again  will  be  nearly  at  atmospheric  pressure. 

For  the  room  air  sampling  experiments  a  tygon  tube  9.5  mm  (3/8") 
i.d.  by  1  metre  long  led  from  the  classifier  aerosol  inlet  and  opened 
horizontally  in  the  room  at  a  height  of  seven  feet.  No  further 
precautions,  e.g.  to  secure  isokinetic  sampling,  wore  felt  necessary  for 
this  preliminary  work. 


ELECTROSTATIC  CLASSIFIER 


The  theory  of  this  Instrument  has  been  developed  by  Hoppel  and 
others  (8,9)  and  the  expressions  relating  to  its  use  as  a  mobility  analyzer 
collected  In  convenient  form  in  a  TSI  publication  by  Keady  (7).  The 
descriptions  and  derivations  given  in  what  follows  are  drawn  entirely  from 
these  sources,  and  from  the  operating  manual  for  the  classifier.  In  brief. 
It  has  the  form  of  a  vertical  cylindrical  condenser,  with  inner  and  outer 
concentric  stainless  steel  cylinders  between  which  the  aerosol  passes. 

Clean  (sheath)  air  is  Introduced  at  the  top,  and  caused  tc  flow  laminarly 
down  the  space  between  outer  cylinder  and  hollow  central  rod.  The  aerosol 
sample  stream  is  Introduced  as  a  thin  annular  layer  along  the  outer  wall, 
avoiding  turbulent  mixing  with  the  sheath  air.  Under  the  influence  of  an 
electric  field  applied  between  outer  (*ve)  and  inner  (**ve)  cylinder. 
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positively  charged  particles  are  attracted  toward  the  centre  rod  as  the  air 
streams  downward.  A  narrow  opening  around  the  perimeter  of  the  central 
rod,  near  its  base,  permits  withdrawal  of  the  fraction  of  aerosol  of  such 
mobility  that,  with  voltage  and  air  flows  set,  and  the  fixed  dimensions  of 
the  apparatus,  it  just  reaches  the  central  rod  at  this  point.  The  excess 
air  is  removed  by  a  vacuum  pump  which  maintains  air  flow  in  the  system. 
Figure  1  shows  the  flows  with  the  terminology  used  here.  Qa  is  evidently 
only  a  small  fraction  of  the  total  air  flow,  Qt.  It  is  also  apparent  that 
the  flows  are  connected  by  the  equation: 


Qe  ♦  «a 


[1] 


Qc,  Qe  and  Qa  are  controlled  by  valves  and  these  flows  are  measured  by  mass 
flowmeters.  It  is  also  necessary,  in  order  to  set  air  flows  to  their 
desired  values,  to  install  a  valve  on  the  sample  inlet  to  control  Q9  as 
well.  To  make  calculations  possible,  the  transfer  function  (to  be  defined 
below)  must  be  symmetrical  and  triangular  in  form,  as  shown  in  Figure  2(a), 
and  for  this  Qa  must  be  set  equal  to  Q3,  from  which  it  follows  that  Qq  - 
Qe.  As  in  fact  the  two  latter  flows  are  considerably  larger  than  the 
former,  the  rough  adjustment  is  made  in  reverse  order,  as  follows. 

With  the  clean  air  inlet  opened  wide,  the  vacuum  pump  is  turned  on 
and  the  excess  air  valve  opened  or  closed  until  excess  air  flow  is  at  the 
desired  rate.  Inlet  (clean)  air  flow  is  now  approximately  the  same  as 
excess  air  flow.  The  CNC  pump  is  turned  on  and  the  monodisperse  aerosol 
valve  is  now  opened  to  a  flow  of  0.25-0.30  1/rain,  approximately  the  intake 
of  the  CNC  (0.3  1/raln).  The  sample  (polydisperse)  valve  is  now  opened  to 
provide  an  adequate  sample  flow,  and  controls  are  then  triraraed  all  round 
until  Qa  «*  Q  (at  the  desired  value,  in  the  range  2.5  -  5  1/mln)  and 
accordingly  Qa  -  Q3  at  the  considerably  lower  value  '0.25  -  0.30  1/min. 

The  larger  the  particle,  the  lower  its  electrical  mobility,  and  the 
longer  it  will  require  to  migrate  from  outer  to  inner  cylinder.  On  the 
other  hand,  the  lower  the  main  flow  rates  Qe  and  Qc,  the  more  time  is 
available  for  this  transport.  Ideally  these  flow  rates  are  set  at  values 
that  will  permit  all  the  aerosol  up  to  the  largest  particles  to  be  detected 
as  the  voltage  is  varied  from  aero  to  its  maximum  value,  near  11,000 
volts. 


It  is  clearly  necessary  that  the  charge  distribution  of  the  aerosol 
be  known,  so  that  for  a  given  sice,  when  the  concentration  in  the  *1  charge 
state  is  determined,  the  concentration  of  all  particles  in  all  charge 
states  be  calculable. 


3! 


3? 


2, 


<b> 


Figure  2:  (a)  Transfer  function  for  classifier 

conditions  described  in  text. 

(b)  Depiction  of  contiguous  transfer 
with  proper  choice  of  classifier 


under 

functions, 

voltage. 
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A  Xr**  neutralizer  in  the  inlet  sample  stream  reduces  the  original 
random  charge  distribution  of  the  aerosol  to  that  corresponding  to  the 
Boltzmann  equilibrium. 

A  recent  calculation  (10)  of  this  distribution  as  a  function  of 
particle  size  appears  in  Table  I.  As  expected,  the  distribution  is 
symmetrical,  with  the  uncharged  component  at  each  size  always  being  the 
largest  single  fraction.  By  plotting  these  results  as  in  Figure  3  one  can 
interpolate  to  get  the  distribution  at  other  diameters.  It  is  to  be  noted 
that  because  of  the  polarity  of  the  classifier,  only  positively  charged 
particles  are  selected;  also,  that  for  a  heterodisperse  aerosol  it  is 
assumed  that  the  particles  in  each  si2e  range  adopt  the  distribution 
appropriate  to  that  si2e,  all  these  distributions  coexisting  in  the 
aerosol . 


TABLE  I 

Distribution  of  Charges  on  Aerosol  Particles 
According  to  Boltzmann's  Law 


D 

tun) 

Percent  of  particles  carrying 

n  elementary  charge  units 

np  * 

-9 

-3 

-2 

-1 

0 

♦  1 

♦  2 

♦3 

♦  9 

0.01 

0.3H 

99.32 

0.3H 

0.02 

5.23 

89.53 

5.23 

0.09 

0.23 

16.22 

67.10 

16.22 

0.23 

0.06 

0.01 

1 .25 

21.30 

59.  SS 

21.30 

1.25 

0.01 

o.oe 

0.08 

2.78 

23-3' 

97.53 

23.37 

2.78 

0.08 

0.10 

0.26 

9.39 

29.09 

92.52 

29.09 

9.39 

0.26 

0.20 

0. 

32 

2.33 

9.66 

22.63 

30.06 

22.63 

9.66 

2.33 

0.32 

0.90 

2. 

19 

5.92 

12.05 

13.  99 

21.26 

16.M 

12.05 

5.92 

2.19 

0.60 

3. 

82 

7. HI 

11.89 

15.79 

17.36 

15.79 

11.89 

7.  HI 

3.82 

0.80 

a. 

83 

7.99 

11.32 

19.00 

15.03 

Mi.  00 

n.32 

7.99 

9.83 

1 .00 

5. 

92 

8.06 

10.71 

12.70 

13.  H5 

12.70 

10.71 

8.06 

5.92 
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Figure  3:  Baltznsann  equUibriua  distribution  of  charge  in 
aerosols,  as  a  function  of  size. 
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In  the  DMPS  arrangement  described  by  TSI,  the  CNC  outlet  leads  back 
into  the  flow  line  through  a  regulating  valve,  so  that  the  CNC  is  part 
of  a  closed  loop.  However,  for  the  CNC  to  operate  properly,  air  flow 
through  it,  monitored  by  an  analog  meter,  must  be  maintained  accurately  at 
5  ml/s  or  0,3  1/min.  With  the  TSI  arrangement  it  was  found  impossible  to 
work  out  any  valve  settings  in  which  the  much  larger  flow  Qe  did  not 
over-ride  the  relatively  small  capacity  pump  that  moves  air  through  the 
CNC.  Accordingly  the  arrangement  of  Figure  1,  with  the  CNC  exhausting  to 
the  atmosphere,  was  used. 


THEORY  OF  CLASSIFIER  OPERATION 


The  electrical  mobility  Z  of  a  particle  in  the  classifier  field  is 
expressed  in  units  of  cm/sec  per  volt/cm  or  velocity  per  unit  field 
strength.  It  has  been  shown  (ref  7)  that  the  mobility  selected  by  the 
classifier  under  a  given  set  of  operating  conditions  is  related  to  the 
other  parameters  by, 


z  ® 


[«t  •  14  <Q3 


+  Qa) j  In 


Hi 

!"S 


2ttVL 
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in  v.'hich  the  flows  are  (in  this  derivation)  expressed  in  cm V sec,  rj  and  ra 
are  the  inner  and  outer  radii  of  the  classifier,  V  the  voltage  and  L  the 
passage  length  from  sample  inlet  to  monodisperse  aerosol  outlet.  Under 
given  conditions,  Z7-K,  and  as  V  is  increased  particles  of  successively 
smaller  mobility  (i.e.  larger  size)  are  selected,  Mobility,  or 
velocity/unit  force  is  expressible  in  another  way  from  the  Stokes  equation 
F  -  3IT1Q.Y,  the  symbols  having  the  usual  significance,  (note  that  small  v 
C 

here  is  velocity,  to  be  distinguished  from  voltage  V).  C  is  the  Cunningham 
correction  which  becomes  important  in  the  cubmicron  region.  -  varies 

F 

inversely  with  diameter;  from  these  equations  a  relation  between  voltage 
and  size  selected  is  at  once  derivable. 

The  selected  mobility  has  a  finite  width 


(®s  +  ®a^  ln  r, 


AZ  * 


2irVL 


[3] 


10  - 


Within  this  width  the  probability  of  selection  of  a  particle  is  expressed 
by  a  transfer  function  which  varies  from  1  at  Z,  the  central  value  to  zero 
at  Z  ±  *r  ,  the  extremes.  It  is  thus  representable  as  a  triangular  area 
(Figure  2(a)).  It  is  also  seen  that  only  half  the  particles  in  the 
mobility  band  AZ  are  selected.  If  this  analysis  is  correct,  then  by 
choosing  voltages  properly  one  can  divide  up  the  aerosol  into  a  series  of 
contiguous  mobility  bands,  and  the  total  measured  particle  count  will  be 
half  the  actual  concentration  of  the  aerosol. 

The  condition  for  contiguous  mobility  bands  is,  referring  to  Figure 
2(b),  that. 


Z2  -  Zi  »  'h  ( AZa  +  AZX) 


Referring  to  auations  2  and  3, 


Z  -  &  where  Ka 
V 


[Qt  -  fc(Qa  ♦  Qa>]  m£i 


2«L 


henoe  Z.  »  ^  ,  Z»  ■ 

V»  Vs 


Kb 

and  AZ  ■  —  where  Kb  • 

V 


(Qs  *  V  ln  FT 

2kL 


and  AZ,  -  ^  ,  AZS  -  ££ 

V,  Va 


Substituting  in  equation  (iJ),  one  finds 


V, 


Ka  -  Kb/2 
Ka  ♦  Kb/2 


which  reduces  to 


V, 


Q  ♦  Q 
s  a 


[4] 


1 


[5) 


a  dimensionless  ratio.  Hence  for  any  set  of  flows  defined  by  Q  ,  Qa  and 
Qt,  a  series  of  voltages  in  geometric  progression  will  survey  the  entire 
aerosol  in  mobility  bands  touching  at  the  base  as  in  Figure  2(b). 


GENERAL  PROCEDURE 


Operation  of  the  classifier  and  CNC  is  described  in  company  manuals. 
The  system  shown  schematically  in  Figure  1  was  set  up,  using  a  Cast  pump  to 
draw  air.  When  measuring  instruments  and  air  flows  had  stabilized,  the 
classifier  was  set  at  a  series  of  voltages  and  the  resulting  CNC  readings 
noted.  This  entailed  recording  by  hand  and  averaging  readings  which 
fluctuated  about  a  final  steady  value.  The  procedure  for  converting  raw 
data  into  finished  results  is  rather  complex;  accordingly  one  example  is 
worked  out  in  nearly  complete  detail  in  Appendix  A.  In  this  example,  the 
TSI  3076  atomizer  was  used  as  source,  with  Q0  “  Qe  -  3  1/min,  and  Qa  »  Qg» 
0.25  1/min,  so  that  Qt  -  3.25  1/min. 

From  equation  5  it  is  seen  that  to  cover  the  entire  aerosol  with 
contiguous  mobility  bands  the  ratio  of  successive  voltages  R  *  1  -  Qg  +  Qa, 

which  In  this  oase  equals  0. 81*61 .  Thus  if  the  final  highest  voltage  is  set 
at  10000V,  the  sequence  will  be  10000,  10000R,  10000R*  etc.  or  10000,  8401 , 
7159,  6057  etc.  volts,  as  shown  in  Table  A-1.  Usually  measurements  were 
taken  with  increasing  voltages,  starting  at  a  value  at  which  count  readings 
first  beoame  significant.  Repeat  runs  and  prolonged  runs  at  a  single 
voltage  were  made  at  intervals,  as  a  check  on  system  stability. 

The  accuracy  of  sizing  was  investigated,  using  highly  monodisperse 
aerosols  of  polystyrene  latex  spheres  (PSL).  These  were  nebulized  in  the 
FT3  *100  and  after  drying  passed  into  the  classifier,  the  voltage  being 
adjusted  until  the  CNC  count  reached  a  maximum.  From  equation  2  the 
mobility  2  was  calculated,  and  hence  the  particle  diameter,  as  explained  in 
Appendix  A.  Measured  diameters  at  four  flow  rates  are  listed  below. 


TABLE  II 


Flow  Rates 

Measured 

Diameter 

2.5  1/min 

i. 00 

.59 

.  42 

.32 

3.0  1/rain 

- 

.58 

.i)2 

.32 

1/raln 

- 

.57 

.U1 

.31 

5.0  1/min 

- 

- 

41 

Nominal  diameter,  ym 

.916 

.562 

.399 

.305 

-  12  - 


Because  of  the  lower  size  cutoff  limits  at  higher  flow  rates,  the 
larger  sizes  become  inaccessible  to  measurement.  Agreement,  apart  from  the 
.916  my  diameter  at  2.5  1/min,  appears  to  be  satisfactory.  The 
concentration  maxima  observed  are  not  as  sharp  as  might  be  expected  from 
the  standard  deviations  on  latex  diameter  quoted  by  the  manufacturer; 
however  these  numbers  are  obtained  by  electron  microscopy,  while  the  PSL 
employed  here  may  have  undergone  some  swelling,  and  is  not  completely  dried 
out  during  aerosol  preparation. 


DYNATECH  OPERATION 


The  3076  atomizer  and  Q127  were  sampled  as  described  earlier.  With 
the  bynatech,  sampling  was  done  through  a  9.5  mm  (3/8")  i.d.  tygon  tube, 
about  3  metres  ir  length,  leading  from  the  centre  of  the  chamber  through 
the  col*  room  to  the  classifier,  which  was  placed  immediately  next  to  the 
cold  room  to  minimize  air  paths. 

In  the  operation  of  this  test  system  air  is  circulated  in  a  closed 
cycle  between  the  control  module  and  the  test  chamber,  the  aerosol  being 
injected  into  the  main  airstream  in  the  test  chamber  inlet  line.  The 
emergen^  air  stream  is  filtered  to  remove  aerosol  prior  to  regeneration,  a 
part  of  the  return  air  being  diverged  through  the  atomizer  and  then 
rejoining  the  main  stream.  '1th  the  cold  room  at  -20°C,  the  return  air  had 
no  opportunity  *-o  warm  up,  and  it  vis  felt  that  this  might  affeot  operation 
of  the  atomizer,  perhaps  by  cooling  the  corn  oil  reservoir  and  affecting 
its  viscosity.  Accordingly  runs  at  -20°J  were  conducted,  first  as 
described  above,  and  then  with  the  return  air  line  in  a  vessel  of  lukewarm 
water' , 


TEST  CONDITIONS 


The  following  runs  were  conducted: 

U) 

Q127  3, 

ii.  5. 

1/min 

(ii) 

TSI  3076 

3.  5 

1/min 

(iii) 

Dynatech 

RT  M9°C) 
O'C 
-20°C 

2.5  3  U  1/min 

3  1/min 

3  1/min  -  no  reheat  of  return  air 

3  1/min  -  with  reheat  of  return  air 


(iv)  Room  air 


2.5  1/min  1530-1630  Weekday 
2030-2130  Saturday 
0830-0930  Weekday 


RESULTS  AND  DISCUSSION 


In  Appendix  A  the  experimental  results  for  TSI  3076  at  a  classifier 
airflow  of  3  1/rain  were  worked  out  in  detail  and  shown  in  Figure  i). 

Results  at  5  1/min  are  similarly  shown  in  Figure  5  and  calculations  for  all 
the  3076  atomizer  work  summarized  in  Table  III (a) .  Comparing  the  figures, 
it  is  seen  that  the  lower  flow  rate,  with  larger  upper  cut-off  size,  covers 
the  aerosol  more  completely.  This  is  reflected  in  the  slightly  larger 
values  for  mass  mean  and  median  diameters.  Second  and  third  order 
corrections  are  also  seen  to  be  smaller  at  3  1/min  (see  Figure  U). 

This  is  presumably  because  at  the  lower  flow  rate  the  upper  size 
cutoff  is  larger  (.8U0  ym  cf  .S^  urn)  and  so  the  assumptions  made  in  the 
calculations  reflect  more  closely  the  actual  situation.  From  Figures  ^  and 
5  it  is  seen  that  the  ratio  concentration  at  largest  size  is,  in  the  first 

peak  concentration 

approximation,  about  5*  at  5  1/min,  but  less  than  0.5<  at  3  1/min. 

Results  for  the  Q127  are  summarized  in  Table  111(b).  Mean  and 
median  diameters  all  lie  in  the  range  0.32  -  0.39  um,  somewhat  larger  than 
the  0.3  um  stipulated.  As  with  the  TSI  3076,  which  also  has  a  reasonably 
monodisperse  aerosol,  the  first  approximation  does  not  differ  much  from  the 
final  corrected  distribution  at  3  1/min.  Results  are  plotted  in  Figure  6. 
It  was  not  felt  worthwhile  to  make  higher  order  corrections  for  the  two 
larger  flow  rates,  4  and  5  1/rain. 

Data  for  the  Dynateoh  appear  in  Table  III(c),  and  the  more 
interesting  results  are  plotted  in  Figures  7  and  8.  Figure  7  shows  the 
first  order  approximation  aerosol  distribution  at  three  temperatures,  the 
experiment  at  -20°C  being  conducted  both  with  and  without  reheat  of  the 
return  air.  Qualitatively  the  shapes  are  similar,  and  the  quantitative 
count  difference  is  evidently  an  artifact  of  the  return  air  temperature. 
This  similarity  is  confirmed  by  examination  of  mean  and  median  diameters 
listed  in  Table  III(c);  in  particular  it  is  seen  on  comparing  results  for  3 
1/rain  at  20#C  and  at  -20°C  with  reheat.  Plots  of  concentration  (1st 
approximation  and  3rd  order  correction)  against  size  for  these  two  runs  are 
nearly  superimposable. 


TABLE  III 


Summary  of  derived  aerosol  properties1 

In  the  first  column,  A  signifies  first  order  approximation, 

B  -  third  order  correction 


Order 

Maximum 

Size 

Measured 

ym 

Maximum 
(Mode) 
in  Number 
Dist  ym 

Number 

Cone 

At  Max 

Number 

Mean 

Diameter 

ym 

Mass 

Median 

Diameter 

ym 

Mass 

Mean 

Diameter 

ym 

Calculated 

Cone 

yg/1 

- 1 - - - 1 

_ _ is] 

1  TSI  3076, Atomizer, 

,  3  1/min 

IMS 

| 

(b)  Q127  3  1/min 


.21 

.23 


Q1 27 .  U  1/min. 


20+5  .22 


Q125.  5  1/min. 


.20 

9.20+5 

.22 

.3* 

\WBM 

.5 

.22-. 2 

■HBfifiil 

.2 

5 

.3 

1 

\wEm 

.3? 

2.5  1/min  R.T. . 


.5*-.55 

.59 


U  1/min  R.T.  . 


.29 

6.58+ 

4 

.29 

3  1/min  0°C  , 


.29 


1/min  -20(>C  No  Reheat 


.26 

5.39+H 

.29 

.50 

.5« 

.52 

.55 

.55 

.59 

Note  1:  See  footnote,  Table  A-1  (Appendix  A),  for  method  of  expressing 
number  concentrations  as  powers  of  10. 
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:  Size  distribution,  TSI  3076  atomizer,  3  1/min, 
several  orders  of  approximation. 
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Figure  5 


Sl2e  distribution,  TSI  3076  atomiser,  5  1/min, 
several  orders  of  approximation. 
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Figure  6 


Size  distribution,  Q127,  3  1/rain,  several  orders  of 
approximation. 


-  18  - 


The  stated  mass  mean  diameter  for  the  Dynatech  is  0.5  -  0.7  pm.  All 
values  listed  in  Table  III(c)  lie  in  the  range  0.50  -  0.63.  This  diameter 
is  considerably  greater  than  the  two  previous  ones,  and  even  at  the  lowest 
practical  flow  rate,  2.5  1/min,  a  considerable  fraction  is  cut  off  at  the 
upper  end  (see  Figure  8).  This  is  reflected  in  a  sizable  difference  in 
calculated  mean  and  median  diameters  on  passing  from  2.5  to  3.0  1/min  and 
the  effect  is  even  more  noticeable  at  H  1/min.  At  room  temperature,  and 
also  at  -20°C  when  the  return  air  is  reheated,  the  calculated  mass 
concentrations  all  lie  in  the  range  25  ±  5  pg/1,  as  stipulated. 

Results  of  the  three  ambient  air  sampling  experiments  are  plotted  on 
Figure  9,  as  first  order  approximations  (against  log  diameter  in  this  case, 
because  of  the  spread  in  size).  The  appearance  of  a  fraction  apparently 
peaking  somewhere  above  1  pm  requires  further  study. 

The  reliability  of  the  mass  concentrations  in  the  last  column  of 
Table  III  is  questionable.  In  principle  this  measurement  should  be 
independent  of  flow  rate,  so  long  as  this  latter  is  not  so  great  that  an 
appreciable  fraction  of  aerosol  at  the  upper  size  end  is  missed.  Results 
for  the  Dynatech,  when  cooling  problems  were  absent,  were  as  expected,  and 
were  independent  of  flow;  however  for  the  TSI  3076  and  the  Q127,  results 
were  sensitive  to  flow  rates.  With  the  Q127,  calculated  concentrations 
varied  from  -  200  pg/1  to  -  50  pg/1  on  passing  from  3  1/min  to  5  1/min,  as 
compared  with  the  stipulated  100  ±  20  pg/1.  A  similar  but  smaller 
variation  was  observed  with  the  TSI  3076.  Using  this  latter  souroe,  a 
direct  determination  of  aerosol  concentration  was  carried  out,  passing  the 
aerosol  stream  at  known  flow  rates  through  a  filter  and  measuring  the 
inorease  in  weight.  At  3  and  5  1/min  the  concentrations  found  by  weight 
inorease  were  80-85  and  52-57  pg/1  respectively,  as  compared  with  110  and 
70  pg/1  reported  in  Table  III.  Perhaps  in  view  of  the  highly  indirect 
method  of  measurement,  effectively  integrating  over  a  series  of  diameter 
intervals,  agreement  within  a  factor  of  two  is  all  that  can  be  expeoted. 

It  is  of  interest  that  summing  the  particle  counts  in  the  atmospheric 
experiments  leads  to  a  value  slightly  larger  than  lOVom*,  about  50-60J  of 
the  count  obtained  when  the  CNC  sampled  the  laboratory  air  directly. 

Figure  9  indicates  that  a  sizable  fraction  i  1  pm  is  being  missed. 


DIAMETER  (*m) 


Figure  7 


Si2e  distribution,  Dynatech,  B  1/min,  at  several 
temperatures. 
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SUMMARY  AND  CONCLUSION 


Use  of  the  electrostatic  classifier  and  nucleus  counter  in  tandem,  a 
well-established  technique  for  determining  aerosol  size  distributions, 
provided  information  on  the  output  of  three  laboratory  generators  that 
agreed  reasonably  well  with  manufacturer's  claims.  Mean  diameters  for  the 
Q127  and  Dynatech  aerosols  were  close  to  the  values  specified.  Height 
concentrations  were  at  worst  within  a  factor  of  two  of  the  expected  values. 
A  preliminary  examination  was  made  of  an  uncontrolled  laboratory 
atmosphere.  Caution  must  be  exercised  in  interpreting  results  with  broadly 
dispersed  aerosols,  as  from  the  Dynatech.  The  method  described  here  for 
determining  higher  order  corrections  is  convenient  for  occasional  or 
one-time  use. 
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APPENDIX  A 


DETAILS  OF  CALCULATION 


GENERAL 


It  13  desirable  that  the  entire  aerosol  size  range  be  spanned,  so 
that  at  the  highest  classifier  voltage  (corresponding  to  the  lowest 
mobility  range  and  t’.v*  largest  diameter)  the  CNC  count  be  nearly  zero,  or 
at  least  much  smaller  than  the  neak  concentration.  Then  one  is  assured 
that  there  is  on’ y  a  negligible  number  of  still  larger  particles  of  +2  or 
higher  charge  falling  in  these  channels;  i.e.  it  can  be  assumed  that  the 
entire  contents  of  these  channels  are  +1  charged  particles  of  the 
appropriate  .ize.  This  condition  is  approximated  in  the  TSI  3076  and  the 
Q127  atomizers,  and  facilitates  calculation.  As  reen  above,  reducing  the 
air  flow  rate  enables  one  to  push  up  the  upper  size  cutoff  and  record 
larger  peptic] es.  A  practical  limitation  is  the  voltage-flow 
characteristics  of  the  mass  flowmeters;  these  become  increasingly 
sensitive  at  low  flow  rates;  one  cannot  go  below  a  rate  of 
Qc  *  Qe  =•  2.5  1/min,  corresponding  to  a  maximum  diameter  of  -  1.0  ym.  At 
this  flow  rate  the  excess  air  valve  is  barely  cracked  open. 


DETAILED  CALCULATION  FOR  THE  TSI  3076  ATOMIZER,  Qq«  Qe  -  3  1/min 


First  Order  Approximation 


This  assumes  that  only  +1  charged  particles  are  present  throughout, 
and  from  the  raw  CNC  count  at  each  voltage  corresponding  to  a  mobility  and 
therefore  to  a  particle  size,  uses  the  Boltzmann  equilibrium  distribution 
to  calculate  the  total  concentration  at  this  size.  In  the  following  table 


each 

column  has  seven  entries. 

(1) 

Voltage 

(2) 

Raw  CNC  count 

i3) 

Mobility  Z  corresponding  to 
equation  2  above;  here,  for 
0.25  1/min,  Z  -  0.1319'V. 

voltage 
Qc  '  Qe 

(1 ) ,  calculated  from 
-  3  1/min,  and  Qg  «  Qa  « 

A-2 


(4)  Particle  size  (pro)  corresponding  to  mobility.  Values  of  Z  vs  D 
have  been  calculated;  a  table  covering  the  range  D  =  0.001  to 
1.0  ym  appears  in  the  manual  for  the  TSI  3071. 

(5)  The  raw  count  in  (2)  is  doubled,  in  order  to  calculate  the 
total  number  count.  As  pointed  out  above,  the  triangular  form 
of  the  transfer  function  results  in  only  half  the  total  aerosol 
in  each  mobility  channel  being  counted.  In  subsequent 
calculations  of  aerosol  distribution,  mass  mean  diameter,  etc., 
the  number  count  per  channel  derived  from  (5)  appears  twice  in 
compensating  ways,  and  so  these  average  values  are  unaffected 
by  doubling  all  counts. 

(6)  The  %  +1  charge  at  Boltzmann  equilibrium  for  the  si2e  in  (4)  is 
determined  by  interpolating  in  Figure  3. 

(7)  The  raw  count  (5)  is  divided  by  (6)  to  give  a  first 
approximation  to  aerosol  concentration  at  size  (4).  More 
correctly  the  diameter  in  (4)  is  a  diameter  range  corresponding 
to  the  mobility  range  AZ,  and  centred  on  the  nominal  value 
given.  All  plots  in  the  graphs  which  follow,  whether  first 
approximations  or  higher  oraer  corrections,  are  concentrations 
at  a  series  of  diameter  ranges. 

Results  for  the  TSI  3076  are  li3ted  in  the  second  column  of  Table 
A-2  and  plotted  as  the  upper  curve  in  Figure  4.  It  is  apparent  that  to 
this  approximation  the  concentration  at  the  largest  size  measured  00.8  um) 
is  a  small  fraotion  of  the  maximum,  at  0.13  -  0.15  pro. 


Second  Order  Calculation 


The  mobility  of  a  doubly  charged  particle  is  twice  that  of  a  singly 
oharged  particle  of  the  same  size.  As  mentioned  above,  it  is  assumed  of 
necessity  that  the  upper  channels  (diameters)  are  populated  only  by  *\ 
charged  particles  of  the  stated  diameter. 

Inspecting  Table  A-1 ,  It  Is  seen  that  the  mobility  of  .840  um 
particles  0.319  -  5)  on  doubling  to  (2.636  -  5)  lies  between  the  values 
for  .483  urn  (2.574  -  5)  and  .424  ym  (3.042  -  5),  and  closer  to  the  former. 
Then  ♦ 2  charged  particles  of  .840  m  diameter  will  be  disposed  primarily  in 
these  two  channels,  and  mostly  at  .483  ym.  In  the  approximation  used  here, 
all  the  contribution  of  *2  oharged  particles  of  0.840  my  diameter  is 
attributed  to  this  smaller  size.  Obviously  in  computer  calculations  this 
contribution  would  be  spread  over  several  channels,  taking  account  of  the 
fact  that  mobility  widths  are  a  function  of  diameter. 
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TABLE  A-1 

TSI  3076,  3  1/min.  Derivation  of  first  order  approximation,  concentration 
as  a  function  of  size. 


1. 

93 

110 

130 

153 

181 

2. 

3.4  +  3  ,  . 

7.00  +  3 

1.03  +  i4 

2.05  + 

U 

3.12  + 

14 

3. 

1.U18  -  3U) 

1.199  -  3 

1.014  -  3 

8.62  - 

1J 

7.29  - 

14 

4. 

.OHO 

.0438 

.048 

.0523 

.0572 

5. 

6.8  +  3 

1.40  +  U 

2.06  +  14 

14.10  + 

1) 

6.2H  + 

H 

6. 

16.2  m 

17 

18 

20.2 

20.8 

7. 

14.2  +  4(a) 

8.2  ♦  14 

1.1H  +  5 

2.03  ♦ 

5 

3.00  + 

5 

1. 

21 1J 

253 

299 

3514 

HI  8 

2. 

H.  38  +  It 

6.35  +  n 

8.146  +  4 

1.07  ♦  5 

1.52+5 

3. 

6.16  -  it 

5.21  -  1) 

14.441  -  U 

3.72  -  U 

3.16  -  D 

14. 

.0627 

.0688 

.0755 

.0831 

.0912 

5. 

8.76  +  H 

1.27  ♦  5 

1.U92  +  5 

2. ID  +  5 

3.0D  ♦  5 

6. 

21.5 

22.H 

23.0 

23.5 

23.8 

7. 

H.08  +  5 

5.67  ♦  5 

6.449  +  5 

9.10  *  5 

1.28+6 

1. 

D9D 

58D 

690 

815 

96D 

2. 

2.05  ♦ 

5 

2.52  +  5 

2.82  +  5 

2.98  +  5 

2.9D  ♦  5 

3. 

2.6D  - 

D 

2.26  -  D 

1.91  -  D 

1.618  -  D 

1.37  -  D 

D, 

.101 

.111 

.123 

.136 

.150 

5. 

D.10  + 

5 

5.0D  ♦  5 

5.6D  +  5 

5.96  +  5 

5.88  +  5 

6. 

2D.1 

2D.  1 

2D,  0 

23.8 

23.5 

7. 

1.70  * 

6 

2.09  ♦  6 

2.35  +  6 

2.50  ♦  6 

2.50  -  6 

1. 

1139 

13D6 

1591 

1880 

2222 

2. 

2.82  +  5 

2.5D  ♦  5 

2.10  ♦  5 

1.60  ♦  5 

1.08  +  5 

3. 

1.158  -  D 

9.80  -  5 

8.29  -  5 

7.016  -  5 

5.936  -  1 

D. 

.167 

.187 

.207 

.232 

.260 

5. 

5.6D  +  5 

5.08  ♦  5 

D.20  +  5 

3.20  ♦  5 

2.16  ♦  5 

6. 

23.2 

22.9 

22. D 

21.9 

21.3 

7. 

:  .  *13  +  6 

2.22  ♦  6 

1.88+6 

1 .D6  +  6 

1.01  +6 

1. 

2t_v> 

3104 

3669 

4336 

5125 

2. 

8.32  ♦  D 

6.09  ♦  D 

3.92  +  4 

2.55  ♦  D 

1.24  ♦  4 

3. 

5.023  -  5 

D.2D9  -  5 

3.595  -  5 

3.042  -  5 

2.574  -  1 

D. 

.293 

.330 

.373 

.424 

.483 

5. 

1 .664  ♦  5 

1.218  ♦  5 

7.84  ♦  4 

5.10  ♦  4 

2.48  ♦  4 

6. 

20.7 

19.8 

18.9 

18.0 

17.2 

7. 

8.04  ♦  5 

6.15  ♦  5 

4.15  ♦  5 

2.83  ♦  5 

1.44  ♦  5 

TABLE  A-1  (Cont'd) 
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10000 
6.8  +  2 
1.319  -  5 
.840 

1.36  +  3 
13.7 
1.0  +  4 


(a)  In  what  follows,  to  simplify  the  tabular  presentations,  positive  and 
negative  exponents  of  ten  are  indicated  as  shown.  Thus  1.418  -  3  = 
1.418  x  10"*,  4.2  +  4  -  4.2  x  10“. 


TABLE  A- 2 

TSI  3076,  3  1/min.  Concentration  as  a  function  of  size. 


Size,  ym 

1st  Order 

2nd  Order 

3rd  Order 

.840 

1.0  +  4 

1.0  +  4 

1.0  +  4 

.728 

2.05  +  4 

2.05  +  4 

2.05  +  4 

.633 

3.87+4 

3.87  ♦  4 

3.87  +  4 

.552 

8.04  +  4 

8.04  +  4 

8.04  +  4 

.483 

1.44  +  5 

1.38  ♦  5 

1.38+5 

.424 

2.83+5 

2.70  +  5 

2.70  +  5 

.373 

4.15  +  5 

3.90+5 

3.90  +  5 

.330 

6.15  +  5 

5.66  +  5 

5.66  +  5 

.293 

8.04  +  5 

7.25  +  5 

7.15  +  5 

.260 

1.01  +  6 

8.59  +  5 

8.45  +  5 

.232 

.1.46  +  6 

1.25+6 

1.22+6 

.207 

1.83  ♦  6 

1.58  +  6 

1.53+6 

.187 

2.22  +  6 

1.86  ♦  6 

1.78+6 

.167 

2.43  ♦  6 

2.03  ♦  6 

1.94  +6 

.150 

2.50  +  6 

1.94  +6 

1.82  +  6 

.136 

2.50  +  6 

1.85  ♦  6 

1 .81  +6 

.123 

2.35  +  6 

1.63  ♦  6 

1.50+6 

.111 

2.09  ♦  6 

1.39  +  6 

1.24  +  6 

.101 

1.70+6 

1.08+6 

9.25  ♦  5 

.0912 

1.28  +  6 

7.31  +  5 

5.80+5 

.0831 

9.10  +  5 

4.72  ♦  5 

3.15  +  5 

.0755 

6.49  ♦  5 

3.39  +  5 

2.47  ♦  5 

.0683 

5.67  +  5 

3.62  +  5 

1.53+5 

.0627 

4.08  ♦  5 

2.93  ♦  5 

1.59  +  5 

.0572 

3.00  +  5 

2.39  ♦  5 

3.89  ♦  4 

.0523 

2.03  ♦  5 

1.68+5 

3.37  ♦  4 

.048 

1.14  +5 

7.61  ♦  4 

1,89  ♦  4 

.044 

8.2  ♦  4 

5.8  ♦  4 

- 

.040 

4.2  ♦  4 

2.8  ♦  4 

- 

1.  6057  7159  8461 

2.  6.55  +  3  2.98  +3  1.50+3 

3.  2.178  -  5  1.842  -  5  1.559  -  5 

4.  .552  .633  .728 

5.  1.31  +  4  5.96  +  3  3.0  +  3 

6.  16.3  15.4  14.6 

7.  8.04  +  4  3.87  +  4  2.05  +  4 
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The  relationship  between  Z  and  V  (equation  2)  and  the  fact  that 
voltages  increase  in  geometric  ratio  (equation  5)  together  ensure  that, 
just  as  the  first  channel  (0.840  ym)  is  matched  with  the  fifth  (0.483  ym) 
in  this  calculation,  so  subsequent  channels  will  be  similarly  matched. 
Thus,  referring  to  Table  A— 1 ,  the  value  of  Z  for  0.207  ym  (8.29  -  5)  on 
doubling  will  lie  between  z  for  .136  and  .123  ym,  and  again  closer  to  the 
former,  four  places  down  from  .207  ym. 


In  the  calculation  the  data  is  treated  four  channels  at  a  time, 
starting  from  the  largest  size.  The  first  round  of  calculation  is  shown  in 
Table  A-3.  Columns  have  the  following  significance: 

(1)  particle  diameter; 

(2)  first  order  concentration  from  Table  A-2,  assumed  correct; 


(3)  percentage  +2  charged  particles  at  size  1,  from  Figure  3; 

(4)  the  product  of  (2)  and  (3),  is  the  concentration  of  doubly 
charged  particles  to  be  expected  in  the  next  four  channels 
(.483  to  .330  ym  inclusive); 


(5)  the  doubled  raw  count  for  these  next  four  channels,  taken  from 
Table  A-1; 

(6)  (5)  -  (4)  the  corrected  raw  count  for  these  sizes; 

(7)  the  $  +1  charged  particles  for  these  sizes,  again  taken  from 
Figure  3. 

(8)  (6)  ♦  (7)  is  then  the  corrected  second  order  concentrations  for 
the  sizes  .483  to  .330  ym. 


TABLE  A-3 


Second  order  correction,  TSI 

3076,  3  1/min, 

first  round. 

See  text 

significance  of 

column  heading  numerals. 

1. 

2. 

3. 

4. 

.840  ym 

1.0  ♦  4 

11.2$ 

1.12  ♦  3 

.728 

2.05  *  4 

11.55 

2.37  ♦  3 

.633 

3.87  ♦  4 

11.8 

4.57  ♦  3 

.552 

8.04  ♦  4 

12.05 

9.69  ♦  3 

Corrected 

5. 

6. 

7. 

8. 

Value  for 

2.48  +  11 

2.37  ♦  4 

17.2$ 

1.38  *  5 

.483  ym 

5.10  ♦  4 

4.86  ♦  4 

18.0 

2.70  ♦  5 

.424 

7.84  +  u 

7.38  ♦  4 

18.9 

3.90  ♦  5 

.373 

1.22  ♦  5 

1.12  ♦  5 

19.8 

5.66  ♦  5 

.330 
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Setting  up  now  a  second  table  as  above  with  .483  to  .330  um  as  the  four 
entries  in  column  (1),  and  their  second  order  corrected  values  just  found 
in  column  (2),  the  process  is  repeated  in  batches  of  four  until  one  has 
stepped  through  all  sizes  and  reached  the  second  order  values  listed  in 
column  3  of  Table  A-2,  and  shown  in  the  middle  curve  of  Figure  4. 


Third  Order  Calculation 


This  begins  with  the  observation  that  the  mobility  for  .840  pm 
tripled  (3.957  ~  5)  lies  between  mobilities  for  the  seventh  and  eighth 
sizes  down  (.373  and  .333  um)  and  somewhat  closer  to  the  eighth  (ie.  seven 
places  down).  Accordingly  the  data  is  treated  in  batches  of  seven.  The 
first  round  of  calculation  is  shown  in  Table  A-4.  Contributions  of  +3 
charged  particles  from  diameters  .840  to  .373.  in  column  (1),  are 
attributed  to  the  next  seven  diameters  .330  to  .167  um. 


TABLE  A-4 

Third  order  correction,  TSI  3076,  3  1/min,  first  round.  See  text  for 
significance  of  column  heading  numerals. 


1. 

2. 

3. 

4. 

.840 

um 

1.0  +  4 

8. 0H 

8.0  + 

2 

.728 

2.05  ♦  4 

7.8 

1.60  + 

3 

.633 

3.87  *  4 

7.55 

2.92  + 

3 

.552 

8.04  +  4 

7.2 

5.79  + 

3 

.483 

1.38  ♦  5 

6.65 

9.18  ♦ 

3 

.424 

2.70  +  5 

6.1 

1.65  ♦ 

4 

.373 

3.90  ♦  5 

5.5 

2.15  ♦ 

4 

CORRECTED 

5. 

6. 

7. 

8. 

VALUE  FOR 

1.12 

5 

1.12  ♦  5 

19.8* 

5.66  + 

5 

.330  um 

1 .50 

♦ 

5 

1.48  ♦  5 

20.7 

7.15  ♦ 

5 

.293 

1.83 

♦ 

5 

1.80  ♦  5 

21.3 

8.45  ♦ 

5 

.260 

2.73 

♦ 

5 

2.67  ♦  5 

21 .9 

1.22  «• 

6 

.232 

3.53 

♦ 

5 

3.43  ♦  5 

22.4 

1.53  * 

6 

.207 

4.26 

♦ 

5 

4,07  ♦  5 

22.9 

1.78  + 

6 

.187 

4.71 

♦ 

5 

4.28  ♦  5 

23.2 

1.94  ♦ 

6 

.167 

Column  (2)  lists  the  best  available  concentrations  for  (1),  the 
first  approximation  for  the  four  largest  sizes,  and  for  the  remaining 
three,  the  seoond  order  corrected  value  found  above.  (3)  is  now  the 
percentage  +3  charge  for  the  sizes  in  column  (1)  taken  from  Figure  3. 

Column  (5),  the  raw  CNC  counts  are  the  best  available  values,  namely  the 
corrected  raw  counts  from  the  second  order  determination  just  made.  The 
remaining  columns  have  the  same  significance. 

This  is  carried  through  in  batches  of  seven.  As  before,  the  third 
order  corrections  of  column  (8)  are  carried  forward  to  the  second  round  and 
placed  in  column  (2).  Table  A-2  column  and  the  lower  curve  of  Figure  4 
show  the  third  order  corrected  values. 


Observations 


Down  to  a  size  -0.1  pm  (see  Figure  4),  the  third  order  correction  is 
hardly  to  be  distinguished  from  the  second,  and  over  most  of  the  range  >0.1 
pm  both  of  these  are  only  minor  corrections  to  the  first  approximation. 
There  are  3 elf- compensating  features  resulting  from  the  skewed  normal 
distribution  observed  in  Figure  4.  At  the  upper  size  end  the  %  content  of 
multiply  charged  particles  is  larger  but  these  corrections  are  applied  to 
sizes  whose  population  is  considerably  greater  to  begin  with,  and  make 
little  difference.  By  the  same  token,  corrections  become  significant  in 
the  lower  size  range,  on  the  left  hand  side  of  the  ourve,  but  here  the  mass 
(though  not  the  number)  contribution  to  the  total  aerosol  is  small;  in  the 
applications  of  this  instrument  the  mass  of  material  penetrating  a  filter 
is  the  most  important  factor. 

The  method  described  above  is  open  to  the  objection  that,  by 
applying  all  seoond  order,  and  then  all  third  order  corrections,  in 
sequence,  certain  interactions  are  overlooked.  This  was  investigated, 
taking  selected  diameters  and  setting  up  a  tree  in  which  these  are 
connected  to  larger  diameters  by  second  or  third  order  effects,  stepping 
back  four  or  seven  places  respectively,  until  one  roaches  one  or  the  other 
of  the  four  largest  diameters,  whose  concentrations  are  by  definition 
assumed  correct. 

Suoh  a  tree  is  shown  below.  Starting  from  the  four  largest  sizes 
(enclosed  in  boxes)  one  oan  apply  all  second  and  third  order  corrections  to 
the  raw  count  data  for  smaller  sizes. 
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2 


.0755 


The  results  of  this  calculation  are  shown  for  a  number  of  sizes  in  the  last 
column. 


TABLE  A-5 


3rd  ORDER 


SIZE,  ym 

1st  ORDER 

2nd  ORDER 

3rd  ORDER 

REFINED 

.330 

6.15  ♦  5 

5.66  +  5 

5.66  +  5 

5.62  +  5 

.260 

1.01  +6 

8.59  +  5 

8.45  +  5 

8.47  +  5 

.232 

1,46  +  6 

1.25+6 

1.22+6 

1.22+6 

.167 

2.43  ♦  6 

2.03  ♦  6 

1.94  +  6 

1.94  +  6 

.150 

2.50  +  6 

1.94  ♦  6 

1.82+6 

1.85+6 

.111 

2.09  ♦  6 

1.39  +  6 

1.24  +  6 

1.27  +  6 

.0755 

6.49  +  5 

3.39  ♦  5 

2.47  +  5 

2.79  +  5 

It  is  evident  that  the  sequential  method  slightly  over-estimates  the 
3rd  order  correction,  but  in  view  of  the  measurement  precision  the 
difference  is  insignificant,  and  the  simple  third  order  correction  is 
entirely  adequate  even  at  sizes  (.0755  ym)  where  this  correction  is  large, 
and  modifies  considerably  the  first  order  approximation. 
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Calculation  of  Size  Distributions 


From  a  set  of  data,  number  concentration  vs  diameter,  by  which  as 
always  is  understood  aerosol  concentrations  in  a  series  of  diameter 
increments,  are  derived  other  measures  of  aerosol  distribution,  as  is  shown 
in  Table  A-6.  Here  columns  A  and  B  are  diameters  and  third  order  corrected 
concentrations  for  the  TSI  3076,  Derivations  in  the  succeeding  columns  are 
explained  in  the  notes  appended  to  the  table. 


TABLE  A-6 

Calculation  of  aerosol  properties,  TSI  3076,  3  1/min 


M 

B 

D 

E 

F 

G 

Efil 

3rd  Order 

(Size) 1 

B  X  C 

Cumulative 

A  x  D 

A  x  B 

■X 

Cone 

(um) 3 

-  mi 

D  -  M 

'  midi 

-  njdj 

.840 

1.0 

+• 

4 

.59270 

5927 

213434 

4987 

Ml 

.728 

2.05 

4 

.38583 

7910 

207507 

5758 

.633 

3.87 

+ 

4 

.25364 

9816 

199597 

6214 

24497 

.552 

8.04 

4 

.16819 

13522 

189781 

7464 

44381 

.483 

1.38 

•f 

5 

.11268 

15550 

176259 

7510 

66654 

.424 

2.70 

5 

.076225 

20580 

160709 

8727 

114480 

.373 

3.90 

+ 

5 

.051895 

20239 

139859 

7549 

145470 

.330 

5.66 

+ 

5 

.035937 

20340 

11  9620 

6712 

186780 

.293 

7.15 

+ 

5 

.025154 

17985 

99280 

5270 

209495 

.260 

8.45 

5 

.017576 

14852 

81295 

3861 

21  9700 

.232 

1,22 

6 

.012487 

1  5234 

66443 

3524 

283040 

.207 

1.53 

6 

.008869 

13570 

51209 

2809 

316710 

.187 

1.78 

+ 

6 

.006539 

11640 

37639 

2177 

332860 

.167 

1.94 

-f 

6 

.004657 

9035 

25999 

1509 

323980 

.150 

1.82 

♦ 

6 

.003375 

61  42 

16964 

921 

237000 

.136 

1.81 

* 

6 

.002515 

4552 

10822 

619 

2461  60 

.123 

1.50 

+ 

6 

.0018608 

2791 

6270 

343 

184500 

.111 

1.24 

6 

.001 368 

1696 

3479 

188 

137640 

.101 

9.25 

* 

5  ! 

.001030 

953 

1783 

96 

93425 

.0912 

5.80 

5 

.000759 

440 

830 

40 

52896 

.0831 

3.15 

+ 

5  1 

.000574 

181 

390 

15 

26176 

.0755 

2.4? 

+ 

5 

.000430 

106 

209 

8 

18648 

.0688 

1.53 

+ 

5 

.000326 

50 

103 

3 

10526 

.0627 

1.59 

+ 

5 

.000246 

39 

53 

2 

9969 

.0572 

3.89 

+ 

4 

.0001 87 

7 

14 

- 

2225 

.0523 

3.37 

+ 

4 

.000143 

5 

7 

- 

1702 

.048 

1.89 

4- 

4 

.00011 

2 

2 

- 

907 

IB  -  1  8.4  x  10* 


ID  -  -M 


I  m1d1  -  76306 


I  nid1  =  3.349,000 
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Mass  median  diameter  Entries  in  Column  C  are  cubes  of  the  diameters  listed 
in  Column  A,  so  that  quantities  derived  from  this  column  are  proportional 
to  volume  and  therefore  to  mass.  Entries  in  Column  D,  the  product  of  the 
corresponding  lines  in  columns  B  and  C,  are  then  proportional  to  the  mass 
contained  in  each  increment  of  diameter.  The  next  column  lists  the 
cumulative  value  of  column  D,  beginning  the  addition  at  the  bottom  of  the 
column,  so  that  the  first  entry  (21  3,  M3  4)  is  proportional  to  total  aerosol 
mass.  From  this,  the  mass  median  diameter  can  be  found,  the  diameter  above 
and  below  which  aerosol  masses  are  equal.  213,434  *  2  •  106,717  and  the 
size  corresponding  to  this  halfway  point  is  0.30  pm. 

Mass  mean  diameter  is  defined  as  midi.  Column  F  is  obtained  by 

multiplying  corresponding  lines  of  A  (dj)  and  D  (m,),  The  sum  of 
column  is  73.606,  which  on  division  by  the  cumulative  value  of  E, 
gives  0.36  pm  as  the  mass  mean  diameter. 

Similarly  the  number  mean  diameter  is  defined  as  Enjdj  in 

column  G,  obtained  by  multiplying  columns  A  and  B  gives  values  of 
Its  sum  is  divided  by  total  number  count,  the  sum  of  column  B,  to 
um. 

Lastly,  assuming  the  derivation  of  equation  5  is  oorrect, 
the  mobility  bands  are  contiguous,  an  estimate  of  the  total  mass 
concentration  (wt/vol)  of  the  aerosol  can  be  made.  The  sum  of  oolumn  D  is 
proportional  to  the  total  aerosol  mass.  Inserting  the  proper 
proportionality  factor  to  convert  to  units  of  pg/1,  and  the  density  of  D0P 
to  convert  volume  to  mass,  one  finds  213,434  x  g  x  10"*  x  0.98  •  110  ug/1. 

Remarks  on  Significant  figures  The  diameters  of  column  A  are  carried  to 
three  significant  figures,  one  more  than  justified,  mainly  to  facilitate 
the  juxtaposing  of  channels  when  calculations  involving  2Z  and  3Z,  for 
higher  order  corrections,  are  carried  out.  The  concentrations  of  column  B 
are  good  at  best  to  two  figures. 

The  purpose  in  carrying  the  extra  numbers  in  column  D  to  0  is  to 
show  the  contrast  between  D-F  and  G.  The  first  three  of  these,  depending 
on  volume,  rapidly  become  insignificant  below  a  diameter  of  0.1  pm,  so  that 
these  entries  could  in  fact  be  dropped  without  affeoting  the  results. 

Column  G  depends  on  number,  not  mass,  and  the  contribution  from  particles 
below  o.l  pm  is  still  appreciable.  The  derived  aerosol  measures  are  all 
rounded  off  to  two  figures. 

It  is  evident  that  for  purposes  of  comparison  one  can  often  use  the 
first  order  approximation  or  even  the  raw  data. 
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